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Abstract 



The use of the differential gluon structure function of the proton J-{x, Q^) introduced 
. by Fadin, Kuraev and Lipatov in 1975 is called upon in many applications of small-x 

^1 I QCD. We report here the first determination of J-{x,Q'^) from the experimental data 

' on the small-x proton structure function F2p{x,Q'^). We give convenient parameteriza- 

O ' tions for J^{x,Q'^) based partly on the available DGLAP evolution fits (GRV, CTEQ 

& MRS) to parton distribution functions and on realistic extrapolations into soft re- 



I gion. We discuss an impact of soft gluons on various observables. The x-dependence of 

the so-determined J-{x,Q'^) varies strongly with and does not exhibit simple Regge 
Oh. properties. None the less the hard-to-soft diffusion is found to give rise to a viable 

(— I I approximation of the proton structure function F2p{x,Q^) by the soft and hard Regge 

components with intercepts Asojj = and ^hard ~ 0.4. 

X. 

■ 1 Introduction: Why unintegrated gluon structure func- 
tions? 

The familiar objects from Gribov-Lipatov-Dokshitzer-Altarelli-Parisi (DGLAP) evolution de- 
scription of deep inelastic scattering (DIS) are quark, antiquark and gluon distribution func- 
tions qi{x,Q'^),q{x,Q^), g{x,Q^) (hereafter x,Q^ are the standard DIS variables). At small 
X they describe the integral flux of partons with the lightcone momentum x in units of the 
target momentum and transverse momentum squared < and form the basis of the highly 
sophisticated description of hard scattering processes in terms of collinear partons On 
the other hand, at very small x the object of the Balitskii-Fadin-Kuraev-Lipatov evolution 
equation is the differential gluon structure function (DGSF) of the target 

•^^'^^ 9iogg2 
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(evidently the related unintegrated distributions can be defined also for charged partons). For 
instance, it is precisely DGSF of the target proton which emerges in the familiar color dipole 
picture of inclusive DIS at small x and diffractive DIS into dijets [Q. Another familiar 
example is the QCD calculation of helicity amplitudes of diffractive DIS into continuum 0, Q 
and production of vector mesons 0, DGSF's are custom-tailored for QCD treatment of 
hard processes, when one needs to keep track of the transverse momentum of gluons neglected 
in the standard collinear approximation 0. 

In the past two decades DGLAP phenomenology of DIS has become a big industry and 
several groups — notably GRV [1^, CTEQ |Tl]] & MRS and some other |jl3[ — keep 



continuously incorporating new experimental data and providing the high energy community 
with updates of the parton distribution functions supplemented with the interpolation routines 
facilitating practical applications. On the other hand, there are several pertinent issues — the 
onset of the purely perturbative QCD treatment of DIS and the impact of soft mechanisms 
of photoabsorption on the proton structure function in the region of large being top ones 
on the list — which can not be answered within the DGLAP approach itself because DGLAP 
evolution is obviously hampered at moderate to small Q^. The related issue is to what extent 
the soft mechanisms of photoabsorption can bias the dependence of the proton structure 
function and, consequently, determination of the gluon density from scaling violations. We 
recall here the recent dispute [jl4| on the applicability of the DGLAP analysis at < 2-4 
GeV^ triggered by the so-called Caldwell's plot |jl5[. Arguably the K-factorization formalism of 
DGSF in which the interesting observables are expanded in interactions of gluons of transverse 
momentum n changing from soft to hard is better suited to look into the issue of soft-hard 
interface. At last but not the least, neglecting the transverse momentum k of gluons is a 
questionable approximation in evaluation of production cross sections of jets or hadrons with 
large transverse momentum. It is distressing, then, that convenient parameterizations of 
DGSF are not yet available in the literature. 

In this communication we report a simple phenomenological determination of the DGSF 
of the proton at small x. We analyze x and dependence of the proton structure function 
F2p{x, Q^) in the framework of the K-factorization approach, which is closely related to the 
color dipole factorization. In the formulation of our Ansatz for J^{x, h?) we take advantage 
of large body of the early work on color dipole factorization 0, |T6], and follow a very 
pragmatic strategy first applied in ^j: (i) for hard gluons with large n we make as much use 
as possible of the existing DGLAP parameterizations of G{x, h?), (ii) for the extrapolation of 
hard gluon densities to small we use an Ansatz Q which correctly describes the color gauge 
invariance constraints on radiation of soft perturbative gluons by colour singlet targets, (iii) 
as suggested by color dipole phenomenology, we supplement the density of gluons with small 

by nonperturbative soft component, (iv) as suggested by the soft-hard diffusion inherent to 
the BFKL evolution, we allow for propagation of the predominantly hard-interaction driven 
small-x rise of DGSF into the soft region invoking plausible soft-to-hard interpolations. The 
last two components of DGSF are parameterized following the modern wisdom on the infrared 
freezing of the QCD coupling and short propagation radius of perturbative gluons. Having 
specified the infrared regularization, we can apply the resulting ^{x^ k^) to evaluation of the 
photoabsorption cross section in the whole range of small to hard Q^. 

The practical realization of the above strategy is expounded as follows: The subject of 
section 2 is a pedagogical introduction into the concept of DGSF on an example of Fermi- 
Weizsacker- Williams photons in QED. Taking the electrically neutral positronium as a target. 
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we explain important constraints imposed by gauge invariance on DGSF at small k^. In 
section 3 we present the /^-factorization approach, which is the basis of our analysis of small-x 
DIS in terms of DGSF. We also comment on the connection between the standard DGLAP 
analysis of DIS and K-factorization and property of soft-to-hard and hard-to-soft diffusion 
inherent to K-factorization. In section 4 we formulate our Ansatz for DGSF. The results of 
determination of DGSF from the experimental data on the proton structure function F2p{x, Q^) 
and real photoabsorption cross section are presented in section 5. In section 6 we discuss an 
anatomy of the so-determined DGSF in the momentum space and comment on the interplay 
of soft and hard components in DGSF, integrated gluon SF and proton structure function 
F2p{x, Q"^). In section 7 we focus on effective intercepts of x-dependence, and the systematics 
of their change from DGSF to integrated gluon SF to F2p{x,Q'^), which illustrates nicely 
the gross features of soft-to-hard and hard-to-soft diffusion pertinent to BFKL physics. The 
subject of section 8 is a comparison of integrated gluon distributions from K-factorization 
and conventional DGLAP analysis of the proton structure function. As anticipated, the two 
distributions diverge substantially at very small x and small to moderate Q^. In section 9 
discuss in more detail how different observables — the scaling violations dF2p{x, Q'^)/d\ogQ'^, 
the longitudinal structure function Fl{x, Q"^) and charm structure function F2'^{x.Q'^) — probe 
the DGSF. In section 10 we summarize our major findings. 



2 Differential density of gauge bosons: the QED exam- 
ple 



Figure 1: 

photons 



The Fermi- Weizsdcker- Williams diagram for calculation of the flux of equivalent 



For the pedagogical introduction we recall the celebrated Fermi- Weizsacker- Williams ap- 
proximation in QED, which is the well known precursor of the parton model (for the review 
see [0). Here high energy reactions in the Coulomb field of a charged particle are treated as 
collisions with equivalent transversely polarized photons — partons of the charged particle, 
fig.|l]. The familiar flux of comoving equivalent transverse soft photons carrying a lightcone 
fraction x-y ^ 1 of the momentum of a relativistic particle, let it be the electron, reads 



Or 



K?dH? dx^ 



kIY X^ 



(2) 
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Here k is photon transverse momentum and = rrif.x^ is the photon longitudinal momentum 
in the electron Breit frame. The origin of in the numerator is in the current conservation, 
i.e. gauge invariance. Then the unintegrated photon structure function of the electron is by 
definition 

^ _dG^_ dnj _aem ( V . . 

If the relativistic particle is a positronium, fig. destructive interference of electromagnetic 
fields of the electron and positron must be taken into account. Specifically, for soft photons 
with the wavelength A = ^ ^ ap, where ap is the positronium Bohr radius, the electromagnetic 
fields of an electron and positron cancel each other and the flux of photons vanishes, whereas 
for A ^ flp the flux of photons will be twice that for a single electron. The above properties 
are quantified by the formula 




:F,'{x,,>^') = 2^ -f— V{k), (4) 



where the factor 2 is a number of charged particles in the positronium and corresponds to 
the Feynman diagrams of fig. |^a, 0d. The vertex function V{k) is expressed in terms of the 
two-body formfactor of the positronium, 

V{k) = 1 - F2(k, -k) = 1- (P| exp(iK(r_ - r+))|P) , (5) 

where r_ — r+ is the spatial separation of e"*" and e~ in the positronium. The two-body 
formfactor F2(k, — k) describes the destructive interference of electromagnetic fields of the 
electron and positron and corresponds to the Feynman diagrams of fig. ^c, It vanishes for 
large enough k > ap^, leaving us with V{k) = 2, whereas for soft gluons one has 

V{k) oc K^al (6) 

One can say that the law is driven by electromagnetic gauge invariance, which guarantees 
that long wave photons decouple from the charge neutral system. 

Finally, recall that the derivation of the differential flux of transverse polarized photons 
would equally hold if the photons were massive vector bosons interacting with the conserved 
current, the only change being in the propagator. For instance, for the charge neutral source 
one finds 

Recall that the massive vector fields are Yukawa-Debye screened with the screening radius 

Rc = my^ . (8) 

To the lowest in QED perturbation theory the two exchanged photons in figs.|l|, ^ do not 
interact and we shall often refer to (0) as the Born approximation for the differential vector 
boson structure function. One can regard as a minimal model for soft k behavior of 
differential structure function for Yukawa-Debye screened vector bosons. 
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Figure 2: The Fermi-Weizsdcker-Williams diagrams for calculation of the flux of equivalent 
photons in positronium. 



3 The insight into the differential density of gluons 
3.1 Modeling virtual photoabsorption in QCD 

The quantity which is measured in deep inelastic leptoproduction is the total cross section of 
photoabsorption 7*p — X summed over all hadronic final states X, where /x, = ±1,0 are 
helicities of (T) transverse and (L) longitudinal virtual photons. One usually starts with the 
imaginary part of the amplitude A^^ of forward Compton scattering 7*p l^p' ■, which by 
optical theorem gives the total cross cross section of photoabsorption of virtual photons 

Q') = ^ ImA^^ , (9) 

at^ixbj, Q') = , ^ ImApo , (10) 

where W is the total energy in the 7*p c.m.s., nip is the proton mass, is the virtuality of 
the photon and Xbj = Q"^/ {Q^ + W"^ — rrip) is the Bjorken variable. 

In perturbative QCD (pQCD) one models virtual photoabsorption in terms of the multiple 
production of gluons, quarks and antiquarks (fig. ^). The experimental integration over the 
full phase space of hadronic states X is substituted in the pQCD calculation by integration 
over the whole phase space of QCD partons 



mi'dr;, ^ y: \m^\' n T —d'>^. , (n) 

n JO Xi 
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Figure 3: The pQCD modeling of DIS in terms of multiproduction of parton final states. 



where the integration over the transverse momenta of partons goes over the whole allowed 
region 



< < 1^2 

- » - 4 



x] 



Ax 



(12) 



The core of the so-called DGLAP approximation [|l| is an observation that at finite x the 
dominant contribution to the multiparton production cross sections comes from a tiny part of 
the phase space 



1 > a^l > X2... > Xn-l > Xn > X 

< < kI... < kI_^ < P < q 



(13) 



in which the upper limit of integration over transverse momenta of partons is much smaller 
than the kinematical limit (|12|). At very small x this limitation of the transverse phase space 
becomes much too restrictive and the DGLAP approximation is doomed to failure. 

Hereafter we focus on how lifting the restrictions on the transverse phase space changes 
our understanding of the gluon structure function of the nucleon at very small x, that is, very 
large - . In this kinematical region the gluon density g{x, Q"^) is much higher than the density 
of charged partons q{x,Q^),q{x,Q^). As Fadin, Kuraev and Lipatov p[ have shown, to the 
leading log- (LL^) approximation the dominant contribution to photoabsorption comes in 
this regime from multigluon final states of fig. alternatively, to the LL^ splitting of gluons 
into gluons dominates the splitting of gluons into qq pairs. As a matter of fact, for the purposes 
of the present analysis we do not need the full BFKL dynamics, in the K-factorization only the 
qq loop is treated explicitly to the LL^ approximation. In this regime the Compton scattering 
can be viewed as an interaction of the nucleon with the lightcone qq Fock states of the photon 
via the exchange by gluons, fig. 11, and the Compton scattering amplitude takes the form 



A 



A, 



qq 



(14) 



Here aa is the and q, q hehcity A, A dependent lightcone wave function of the photon 
and the QCD pomeron exchange gg-proton scattering kernel Agg does not depend on, and 



conserves exactly, the g, g helicities, summation over which is understood in (|14D. 
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Figure 4: The K-factorization representation for DIS at small x. 



The resummation of diagrams of fig. |^ defines the unintegrated gluon structure function 
of the target, which is represented in diagrams of fig. ^ as the dashed blob. The calculation 
of the forward Compton scattering amplitudes (A = 0) is straightforward and gives the k- 
factorization formulas for photoabsorption cross sections [p^ pO|] 

TI . ■' JO J J 



f 



\z' + (1 - zf] 



k — K 



+m 



(15) 



(^L{xBj,Q'^) = Ve^. dz d^k / —-as{q'^)J^{xg,K'^ 

TT J, Jo J J 



4Q'z\l-zf 



1 



(16) 



^k2 + £2 (k-K)2+£2y 

Here ruf and e/ are the mass and charge of the quark f = u,d, s, c, b, .., 

= z{l - z)Q^ + mj , (17) 
the QCD running coupling asiq"^) enters the integrand at the largest relevant virtuality, 



q = maxje + k , k } 



and the density of gluons enters at 



_ + _ (..Ml 
+ Q2 ^"^^ y + Q2 



(19) 



Here Mt is the transverse mass of the produced qq pair in the photon-gluon fusion 7*5? — >■ qq\ 



ml + + (k 



i2 



m; = + — ^ . (20) 

1- z z 

For longitudinal photons only transitions 7^ q\q\ into states with A + A = are allowed. 
In (Tt the terms oc are the contribution of states with A + A = yU, whereas the dominant 
contribution in the scaling regime of ^ comes from the transitions 7t — > qxqx into 
states A + A = 0, when the helicity of the photon is transferred to the angular momentum of 
the quark-antiquark pair. The corresponding transition amplitudes are oc k, k ± k, for more 
discussion see [0]. 

No restrictions on the transverse momentum in the qq loop, k, and gluon momentum, k, 
are imposed in the representation ([l5|), (0). This representation was contained essentially 
in the classic Fadin, Kuraev, Lipatov papers of mid-70's, in the recent literature it is 
sometimes referred to as the K-factorization. 

We note that Eqs. ([TBI) , (|l^) are for forward diagonal Compton scattering, but similar 
representation in terms of the unintegrated gluons structure function holds also for the off- 
forward Compton scattering at finite momentum transfer A, off-diagonal Compton scattering 
when the virtualities of the initial and final state photons are different, 7^ Q^, including 
the timelike photons and vector mesons, = — m^, in the final state. 

The photoabsorption cross sections define the dimensionless structure functions 



i^T,L(x6,,Q') = ^^^aT,L (21) 



and F2 = Ft + Fl, which admit the familiar pQCD parton model interpretation 

Frix,,, Q') = ^]M^b,, Q') + Qfix,,, Q^)] , (22) 

f=u,d,s,c,b,.. 

where qf{xbj, Q^), Qfixbj, Q^) are the integrated density of quarks and antiquarks carrying the 
fraction Xbj of the lightcone momentum of the target and with transverse momenta < Q. 
Hereafter we suppress the subscript "bj" . 

3.2 Where At^-factorization meets DGLAP factorization 

Recall the familiar DGLAP equation |l|] for scaling violations at small x, 

dF2{x,Q^) ^ 2«5(Q^) /"^ r r 2 , /I \2i^/^ n^^ ^^iQ^) ^ n2\ \^ 2 /oQ^ 
-dl^ = ^ ^/^^ L ^^[^ + - ^) ' ^ ^ " ^ ^ ^ ' ^23) 

where for the sake of simplicity we only consider light flavours. Upon integration we find 

F2{x,Q )^}_^ej G{2x,Q). (24) 

f ■'^ Q 37r 

In order to see the correspondence between the Kx-factorization and DGLAP factorization it 
is instructive to follow the derivation of (|2^) from the K^-representation ([15|) . 
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First, separate the ^^-integration in (15) into the DGLAP part of the gluon phase space 



< Q = + k"^ and beyond-DGLAP region k,^ >,Q ■ One readily finds 



k^ + s^ (k-^r + ey U-^^ , if ^'>Q' 



(25) 



n2 



Consider first the contribution from the DGLAP part of the phase space ^Q^- Notice 
that because of the factor in (ESl), the straightforward integration of the DGLAP com- 

2 2 

ponent yields G{xg, Q ) and Q is precisely the pQCD hard scale for the gluonic transverse 
momentum scale: 



V g Q Q J 

The contribution from the beyond-DGLAP region of the phase space can be evaluated as 

'^2 ^"S(g ):F{Xg, ) =2 =4 = =4 =6 ^i^9^ Q )HXg, Q ) 

Q 1^ \Q Q J \Q Q J 



'2^2(1 _^)2Q4 2z{l-z)Q^ 1 



+ ^ J ^(xg, g') log C2(xg, q', . (27) 

The latter form of (pT]) allows to conveniently combine (|26D and (^) rescaling the hard scale 
in the GSF 

G{xg,Q') + T{xg,Q')logC2{xg,Q\z) = G{xg,C2{xg,Q\z)Q') . (28) 
Here the exact value of I{xg, Q ) > I depends on the rate of the K^-rise of J^{xg, k,^). At small 

2 2 2 

Xg and small to moderate Q one finds I{xg, Q ) substantially larger than 1 and C2{xg, Q ^z) ^ 
1, see more discussion below in section 9. 



Now change from dk integration to dQ and again split the z^Q integration into the 

■2 



DGLAP part of the phase space Q <^ \Q'^, where either z<^oil — z<^, and the 
beyond-DGLAP region > |g^, where < z < 1. As a result one finds 

. ^2^2|-2A,2i + A4), if Q'^\Q^ 

where 

= dz[z^ + (1 - z)2]2'"(l - ^)™ (30) 



(29) 
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Let C2 be C2{xg, Q , z) at a mean point. Notice also that Mf ~ Q^, so that Xg ~ 2x. Then 
the contribution from the DGLAP phase space of Q can be cast in precisely the form (^) 

F2{x,Q') ^Ve^. /'^'''^^^G(2x,g^). (31) 

^ DGLAP y ^ Jo ?,Tl ^ ' ^ ' 

The beyond-DGLAP region of the phase space gives the extra contribution of the form 

^J2e}^^G{2x,Q'). (32) 

Eqs. ( ^H) and (|32D immediately reveal the phenomenological consequences of lifting the DGLAP 
restrictions in the transverse momenta integration. Indeed, the DGLAP approach respects the 
following strict inequalities 

^2 < and < Q2 (33) 

— 2 

As we just saw, removing the first limitation effectively shifted the upper limit in the Q 
integral to 7^ Q^, while lifting the second constraint led to an additional, purely non- 

DGLAP contribution. Although both of these corrections lack one leading log-Q^ factor they 
are numerically substantial. As a matter of fact, in section 9 we show that C2 ~ 8. 

The above analysis suggests that the DGLAP and K-factorization schemes converge loga- 
rithmically at large Q^. However, in order to reproduce the result (|3lD and (0) for the full 
phase space by the conventional DGLAP contribution ( P^ from the restricted phase space 
(|13D one has to ask for DGLAP gluon density Gpt{x, Q^) larger than the integrated GSF in the 
K-factorization scheme and the difference may be quite substantial in the domain of strong 
scaling violations. 



3.3 The different evolution paths: soft-to-hard diffusion and vice 
versa 

The above discussion of the contributions to the total cross section from the DGLAP and 
non-DGLAP parts of the phase space can conveniently be cast in the form of the Huygens 
principle. To the standard DGLAP leading logQ^ (LLQ^) approximation one only considers 
the contribution from the restricted part of the available transverse phase space (p!3D . The 
familiar Huygens principle for the homogeneous DGLAP LLQ^ evaluation of parton densities 
in the Xbj-Q'^ plane is illustrated in fig. §a: one starts with the boundary condition p{x, Qq) as 
a function of x at fixed Qq, the evolution paths [z, Q^) for the calculation of p{x, Q^) shown 
in fig. 5a are confined to a rectangle x < z <1, Qq < < Q^, the evolution is unidirectional 
in the sense that there is no feedback on the x-dependence of p{x, Qf) from the x-dependence 
of p(a;, Ql) at Ql > Q\. In fig. 5a we show some examples of evolution paths which are 
kinematically allowed but neglected in the DGLAP approximation. Starting with about fiat 
or slowly rising G(x, Qq) one finds that the larger Q^, the steeper the small-x rise of G(x, Q'^). 
At a; <^ 1 the DGLAP contribution from the restricted transverse phase space (|13D no 



longer dominates the multiparton production cross sections, the restriction (|T3D must be lifted 
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-1 X -3 -2 -1 

10 1 ^10 10 10 1 



a) b) 

Figure 5: The Huygens principle for Q^,x evolution of DIS structure functions with (a) 
DGLAP restricted transverse phase space and (b) for the BFKL x evolution without restric- 
tions on the transverse phase space and hard-to-soft & soft-to-hard diffusion. 



and the contribution to the cross section from small and large > can no longer be 
neglected. The Huygens principle for the homogeneous BFKL evolution is illustrated in fig. |^: 
one starts with the boundary condition jF(xo,(5^) ^ function of at fixed xq <^ 1, the 
evolution paths {z, Q'^) for the calculation of p{x, Q"^) are confined to a stripe x < z < xq, in 
contrast the the unidirectional DGLAP evolution one can say that under BFKL evolution the 
small-x behaviour of p{x, Q^) at large is fed partly by the x-dependence of soft p{x, Q"^) at 
larger x and vice versa. The most dramatic consequence of this soft-to-hard and hard-to-soft 
diffusion which can not be eliminated is that at very small x the x-dependence of the gluon 
structure in the soft and hard regions will eventually be the same: 

lim G{x,Q^)=G{Q') (-^'^ . (34) 



.X, 

The rate of such a hard-to-soft diffusion is evidently sensitive to the infrared regularization of 



pQCD, the model estimates show that in the HERA range of x it is very slow |21, 22 . 



4 The Ansatz for differential gluon structure function 

The major insight into parameterization of DGSF comes from early experience with color 
dipole phenomenology of small-x DIS. In color dipole approach, which is closely related to 
K-factorization, the principal quantity is the total cross section of interaction of the qq color 
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dipole r with the proton target 0, ^ ^ 

r) = ^J — ^— p as [max{K , -} j T{x, k ) , (35) 

which for very small color dipoles can be approximated by 

a{x, r) = -—as {^) G (x, — ) , (36) 



where A ^ 10 comes from properties of the Bessel function Jo{z). The phenomenological 
properties of the dipole cross section are well understood, for extraction of a{x, r) from the 
experimental data see 2^. The known dipole size dependence of a{x,r) serves as a 



constrain on the possible K^-dependence of J-'{x, k,"^). 

As we argued in section 3.2, DGLAP fits are likely to overestimate J^hardi^, f^^) at moderate 
H?. Still, approximation ( p6D does a good job when the hardness A/r"^ is very large, and at 
large we can arguably approximate the DGSF by the direct differentiation of available fits 
(GRV, CTEQ, MRS, ...) to the integrated gluon structure function Gpt{x,Q'^): 

^..(x,/.^)-^^ (37) 

O log 

Hereafter the subscript pt serves as a reminder that these gluon distributions were obtained 
from the pQCD evolution analyses of the proton structure function and cross sections of related 
hard processes. 

The available DGLAP fits are only applicable at k,^ > Ql, see table 1 for the values of 
Qc, in the extrapolation to soft region < Ql we are bound to educated guess. To this end 
recall that perturbative gluons are confined and do not propagate to large distances; recent 



fits |2^ to the lattice QCD data suggest Yukawa-Debye screening of perturbative color fields 
with propagation/screening radius Rc ~ 0.27 fm. Incidentally, precisely this value of Rc for 
Yukawa screened colour fields has been used since 1994 in the very successful color dipole 



phenomenology of small-x DIS 0, |T3- Furthermore, important finding of |jT^ is a good 
quantitative description of the rising component of the proton structure function starting 
with the Yukawa-screened perturbative two-gluon exchange as a boundary condition for the 
color dipole BFKL evolution. 

The above suggests that k'^ dependence of perturbative hard J-'hard{x, k^) in the soft region 
< is similar to the Yukawa-screened flux of photons in the positron, cf. eq. (^, with 
aem replaced by the running strong couphng of quarks CFdsii^"^) and with factor Nc instead 
of 2 leptons in the positronium, for the early discussion see 0], 

= CN^^ (^^)'^-(-) ' (38) 

Here fipt = = 0.75 GeV is the inverse Yukawa screening radius and must not be interpreted 
as a gluon mass; more sophisticated forms of screening can well be considered. Following 
[|l^, |l^ |19|, ^ we impose also the infrared freezing of strong coupling: ^^(k^) < 0.82; recently 



the concept of freezing coupling has become very popular, for the review see 

The vertex function Vn^k,) describes the decoupling of soft gluons, /t <^ from color 
neutral proton and has the same structure as in eq. (^. In the nonrelativistic oscillator model 
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for the nucleon one can relate the two-quark form factor of the nucleon to the single-quark 
form factor, 

F2(^,-^) = Fi(^|^..^) . (39) 
To the extent that R"^ <^ the detailed functional form of F2(/5, — k) is not crucial, the 



simple relation (|59|) will be used also for a more realistic dipole approximation 

Fi(k2) = L^. (40) 

The gluon probed radius of the proton and the charge radius of the proton can be somewhat 
different and A ~ 1 GeV must be regarded as a free parameter. Anticipating the forthcoming 
discussion of the diffraction slope in vector meson production we put A = 1 GeV. 

As discussed above, the hard-to-soft diffusion makes the DGSF rising at small-x even in 
the soft region. We model this hard-to-soft diffusion by matching the dependence ( ^8]) to 
the DGLAP fit J-'pt{x, Ql) at the soft-hard interface Ql and assigning to J-'hard{x, k"^) in the 
region of < Ql the same x-dependence as shown by the DGLAP fit J-'pt{x, Ql), i.e.. 



^w(x, ^) = - >^') + ^vtix, '^)e{i^ - Ql) . (41) 



Because the accepted propagation radius Rc ~ 0.3 fm for perturbative gluons is short compared 
to a typical range of strong interaction, the dipole cross section (|35|) evaluated with the DGSF 
( ^Tj ) would miss an interaction strength in the soft region, for large color dipoles. 



In Ref . [[Tq, p!7| interaction of large dipoles has been modeled by the non-perturbative, 
soft mechanism with energy-independent dipole cross section, whose specific form |l^, has 
been driven by early analysis |]19[ of the exchange by two nonperturbative gluons. More 
recently several closely related models for (Jso{t{r) have appeared in the literature, see for 
instance models for dipole-dipole scattering via polarization of non-perturbative QCD vacuum 
[p9|] and the model of soft-hard two-component pomeron In order to reproduce the 

required interaction strength for large dipoles, we introduce the genuinely soft, nonperturbative 
component of DGSF which we parameterize as 



^St(x, k') = asoftCpN,^^ ( ] Vm{^^) . (42) 




The principal point about this non-perturbative component of DGSF is that it must not be 
subjected to pQCD evolution. Thus the arguments about the hard-to-soft diffusion driven 
rise of perturbative DGSF even at small do not apply to the non-perturbative DGSF 
and we take it energy-independent one. Such a nonperturbative component of DGSF would 
have certain high-K^ tail which should not extend too far. The desired suppression of soft 
DGSF at large-K^ and of hard DGSF (^) at moderate and small can be achieved by the 
extrapolation of the form suggested in [|, ^ 

r^) = rl^lix, ^'')-^, + -Fw(x, 1^')-^ (43) 

The above described Ansatz for DGSF must be regarded as a poor man's approximation. The 
separation of small- fc^ DGSF into the genuine nonperturbative component and small-K^ tail 
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of the hard perturbative DGSF is not unique. Specifically, we attributed to the latter the 
same small-x rise as in the DGLAP fits at Ql while one can not exclude that the hard DGSF 
has a small x-independent component. The issue of soft-hard separation must be addressed 
in dynamical models for infrared regularization of perturbative QCD. As we shall see below, 
in section 4.4, the soft component of the above described Ansatz is about twice larger than 
the soft component used in the early color dipole phenomenology |jl6|, 17 1. 



The K-factorization formulas ([T5|) and ([T6| ) correspond to the full-phase space extension 
of the LO DGLAP approach at small x. For this reason our Ansatze for J^hardi^, Q"^) will be 
based on LO DGLAP fits to the gluon structure function of the proton Gpt{x, Q'^). We consider 
the GRV98LO |T0|, CTEQ4L, version 4.6 and MRS LO 1998 pf parameterizations. We 
take the liberty of referring to our Anzatze for DGSF based on those LO DGLAP input as 
D-GRV, D-CTEQ and D-MRS parameterizations, respectively. 



Our formulas ([T^, (H) describe the sea component of the proton structure function. 
Arguably these LL^ formulas are applicable at a; < Xq = 1 ^ 3 • 10^^. At large the 
experimentally attainable values of x are not so small. In order to give a crude idea on finite- 
energy effects at moderately small x, we stretch our fits to x > xq multiplying the above 
Ansatz for DGSF by the purely phenomenological factor (1 — x)^ motivated by the familiar 
large-x behaviour of DGLAP parameterizations of the gluon structure function of the proton. 
We also add to the sea components (|l5), (P^ the contribution from DIS on valence quarks 



borrowing the parameterizations from the respective GRV, CTEQ and MRS fits. The latter 
are only available for > Ql- At x < 10~^ this valence contribution is small and fades away 



rapidly with decreasing x, for instance see ||T7 



5 Determination of the differential gluon structure func- 
tion of the proton 

5.1 The parameters of DGFS for different DGLAP inputs 

Our goal is a determination of small-x DGSF in the whole range of /t^ by adjusting the relevant 
parameters to the experimental data on small-x F2p{x,Q^) in the whole available region of 

as well as the real photoabsorption cross section. The theoretical calculation of these 
observables is based on Eqs. (|15|), (|16|), (^31). 

The parameters which we did not try adjusting but borrowed from early work in the color 
dipole picture are Rc = 0.27 fm, i.e., fipt = 0.75 GeV and the frozen value of the LO QCD 
coupling with Aqcd = 0.2 GeV: 

asiQ') = min 0.82, ^ } . (44) 

(3o log 

K ^^QCD 

We recall that the GRV, MRS and CTEQ fits to GSF start the DGLAP evolution at quite a 
different soft-to-hard interface Qc and diverge quite a lot, especially at moderate and small 
k;^. The value of Qc is borrowed from these fits and is not a free parameter. 

The adjustable parameters are /iso/i, Oso/t, iT^u,d, i^l and n\. We take = m„^rf + 0.15GeV 
and Trie = 1.5 GeV. The role of these parameters is as follows. The quark mass rriu^d defines 
the transverse size of the qq = uu, dd Fock state of the photon, whereas /iJoj^ controls the 



14 



r- dependence of, and in conjunction with a^o/t sets the scale for, the dipole cross section for 
large size qq dipoles in the photon. The both m^^ and fisoft have clear physical meaning 
and we have certain insight into their variation range form the early work on color dipole 
phenomenology of DIS. The magnitude of the dipole cross section at large and moderately 
small dipole size depends also on the soft-to-hard interpolation of DGSF, which is somewhat 
different for different LO DGLAP inputs for Gpt{x, Q^). This difference of DGLAP inputs can 
be corrected for by adjusting nloft ^^^d the value of k\. Because of soft-to-hard and hard-to- 
soft diffusion the DGLAP fits are expected to fail at small we allow for x dependence of 

Evidently, roughly equal values of F2p{x, Q^) can be obtained for somewhat smaller J^{x, Q^) 
at the expense of taking smaller mu,d and vise versa. Therefore, though the quark mass does 
not explicitly enter the parameterization for J-'{x, k^), the preferred value of mu,d turns out 
to be correlated with the Ansatz for DGSF, i.e. each particular parameterization of DGSF 
implies certain mu^d- In what follows we set agoft = 2, = 3.0 GeV^, mu,d = 0.22 GeV, 
so that only k1 and figoft varied from one DGLAP input to another, see Table 1. The soft 
components of the D-GRV and D-CTEQ parameterizations turn out identical. The eye-ball 
fits are sufficient for the purposes of the present exploratory study. The parameters found are 
similar to those used in ^ where the focus has been on the description of diffractive DIS. 



Table 1. The parameters of differential gluon structure function for different DGLAP inputs. 





D-GRV 


D-MRS 


D-CTEQ 


LO DGLAP input 


GRV98LO 10| 


CTEQ4L(v.4.6) Q 


MRS-LO-1998 p| 


Ql GeV^ 


0.895 


1.37 


3.26 


4, GeV2 


1 /2 

(l + 0.0018 log^i) 


(1 + 0.038 log^^)^^^ 


(l + 0.047 log^i)^^^ 


fJ'soft, GeV 


0.1 


0.07 


0.1 



One minor problem encountered in numerical differentiation of all three parameterizations 
for Gpt{x,Q'^) was the seesaw At^-behavior of the resulting DGSF (^7\), which was an artifact 
of the grid interpolation routines. Although this seesaw behavior of DGSF is smoothed out 
in integral observables like G{x, Q"^) or F2p{x, Q^), we still preferred to remove the unphysical 
seesaw cusps and have smooth DGSF. This was achieved by calculation DGSF from (^) at 
the center of each interval of the Q^-grid and further interpolation of the results between 
these points. By integration of the so-smoothed J^pt{x, Q"^) one recovers the input Gpt{x, Q"^). 
The values of Ql cited in Table 1 corresponds to centers of the first bin of the corresponding 
Q^-grid. 



5.2 The description of the proton structure function F2p{x,Q'^) 

We focus on the sea dominated leading log^ region of x < 10~^. The practical calculation 
of the proton structure function involves the two running arguments of DGSF: Xg and k^. 
We recall that in the standard coUinear DGLAP approximation one has -C fc^ <^ and 
Xg ~ 2x, see eq. (p3|). Within the K-factorization one finds that the dominant contribution to 
F2p{x, Q"^) comes from Mf ~ with little contribution from > Q^. Because at small Xg 
the Xg dependence of J-'{xg, Q"^) is rather steep, we take into account the Xg — Xhj relationship 
(p!9[). Anticipating the results on effective intercepts to be reported in section 7, we notice 
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Figure 6: The n-factorization description of the experimental data on F2p{x, Q"^) in the low 
region; black circles are ZEUS BPC data j^^, open triangles denote HI shifted vertex (SV) 
data ^^J, open squares are E665 data j^/ . Solid line represents K-factorization results for 
the D-GRV parameterization of the differential gluon structure function J-'{x, k^). 



that for all practical purposes one can neglect the impact of k on the relationship (|T9D, which 
simplifies greatly the numerical analysis. Indeed, the Xg dependence of T{xg^ k^) is important 
only at large k^, which contribute to F2p{x, Q^) only at large Q^; but the larger Q^, the better 
holds the DGLAP ordering -C /c^, Q^. Although at small to moderate Q"^ the DGLAP the 
ordering breaks down, the Xg dependence of J-'{xg, n^) is negligible weak here. 

As we shall discuss in more detail below, achieving a good agreement from small to moder- 
ate to large is a highly nontrivial task, because strong modification of the soft contribution 
to T{x, Q"^) unavoidably echos in the integral quantity Gr){x, Q"^) throughout the whole range 
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Figure 7: The K-factorization description of the experimental data on F2p{x,Q'^) in the mod- 
erate and high region; black circles and triangles are ZEUS data ^3l\] , ^3^] , open circles 
and triangles show HI data /|^/, j^/ , open squares are E665 data j^/ , stars refer to NMC 
results ^3l{J . Solid line represents n-factorization results for the D-GRV parameterization of 
the differential gluon structure function jF(x, n^). 



of and shall affect the calculated structure function from small to moderate to large Q^. 

The quality of achieved description of the experimental data on the small-x proton structure 
function is illustrated by figs. ^, 0. The data shown include recent HERA data (ZEUS [ pT[| , 



ZEUS shifted vertex H, ZEUS BPC H, HI g, HI shifted vertex ll), FNAL E665 
experiment and CERN NMC experiment [^. When plotting the E665 and NMC data, 
we took the liberty of shifting the data points from the reported values of to the closest 
boxes for which the HERA data are available. For Q"^ < = 0.9 GeV^ the parameterizations 
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Figure 8: A comparison of the n- factorization description of the experimental data on 
F2p{x,Q'^) for several values of with the D-GRV, D-CTEQ and D-MRS parameterizations 
of the differential gluon structure function T{x^ k}\ The contribution to F2p{x, Q"^) from DIS 
off valence quarks is shown separately for larger . 



for valence distributions are not available and our curves show only the sea component of 
F2p(x,(5^), at larger Q"^ the valence component is included. 

At X < 10~^ the accuracy of our D-GRV description of the proton structure function is 
commensurate to that of the accuracy of standard LO GRV fits. In order not to cram the 
figures with nearly overlapping curves, we show the results for D-GRV parameterization. The 
situation with D-CTEQ and D-MRS is very similar, which is seen in fig. where we show 
on a larger scale simultaneously the results from the D-GRV, D-CTEQ and D-MRS DGSFs 
for several selected values of Q^. Here at large we show separately the contribution from 
valence quarks. The difference between the results for F2p{x, Q"^) for different DGLAP inputs 
is marginal for all the practical purposes, see also a comparison of the results for a'^^ for 
different DGLAP inputs in fig. ^ 
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Figure 9: A comparison of the K-factorization description of the experimental data on real 
photoah sorption cross section based on the D-GRV, D-CTEQ and D-MRS parameterizations 
of the differential gluon structure function T{x^h?). The squares show the experimental data 
from 1992-93 direct measurements, the bullets are the results of extrapolation of virtual pho- 
toabsorption to = (j^^ and references therein). The soft component of photoab sorption 
cross section is shown separately. 

5.3 Real photoabsorption cross section a^^. 

In the limiting case of = the relevant observable is the real photoabsorption cross section 
a'^'P. Although the Bjorken variable is meaningless at very small Q^, the gluon variable Xg 
remains well defined at = 0, see eq. (|19|). In fig. |^ we present our results alongside with 



the results of the direct measurements of a"'^ and the results of extrapolation of virtual pho- 
toabsorption cross sections to = 0, for the summary of the experimental data see ||33|. We 



emphasize that we reproduce well the observed magnitude and pattern of the energy depen- 
dence of (t'^p in an approach with the manifestly energy-independent soft contribution to the 
total cross section (which is shown separately in fig. We recall that our parameterizations 
for jF(x,Q^) give identical soft cross sections for the GRV and CTEQ inputs (see table 1). 
The barely visible decrease of fj^j^ towards small is a manifestation of oc (1 — a;)^ large-x 
behaviour of gluon densities. The extension to lower energies requires introduction of the 
secondary reggeon exchanges which goes beyond the subject of this study. 
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In our scenario the energy dependence of a"''^ is entirely due to the x^-dependent hard 
component J-'hard{xg, Q^) and as such this rise of the total cross section for soft reaction can 
be regarded as driven entirely by very substantial hard-to-soft diffusion. Such a scenario has 
repeatedly been discussed earlier [p!6| , |1^, Time and time again we shall see similar effects 
of hard-to-soft diffusion and vise versa. Notice that hard-to-soft diffusion is a straightforward 
consequence of full phase space calculation of partonic cross sections and we do not see any 
possibility for decoupling of hard gluon contribution from the total cross sections of any soft 
interaction, whose generic example is the real photoabsorption. 



6 Properties of differential gluon structure function in 
the momentum space 

6.1 Soft/hard decomposition of DGSF 

Now we focus on the x and behavior of the so-determined DGSF starting for the reference 
with the D-GRV parameterization. The same pattern holds for DGSF based on CTEQ and 
MRS DGLAP inputs, see below. In figs. [1^ and ^ we plot the differential gluon density 



J-'{xg, Q"^), while in fig. |12| we show the integrated gluon density 

Go{x,Q^)= (45) 

JO K 

Here the subscript D is a reminder that the integrated Gr){x,Q^) is derived from DGSF. As 
such it must not be confused with the DGLAP parameterizations Gpt{x, Q"^) supplied with 
the subscript pt. 



Figs. and |TT| illustrate the interplay of the nonperturbative soft component of DGSF 
and perturbative hard contribution supplemented with the above described continuation into 

< various x and h?. The soft and hard contributions are shown by dashed and 

dotted lines respectively; their sum is given by solid line. 

Apart from the large-x suppression factor (1 — x)^ our non-perturbative soft component 
does not depend on x. At a not so small x = 10^^ it dominates the soft region of < 1 ^ 2 
GeV^, the hard component takes over at higher k^. The soft-hard crossover point is close to 
/Xpj but because of the hard-to-soft diffusion it moves with decreasing x to a gradually smaller 



6.2 Soft /hard decomposition of the integrated gluon structure func- 
tion 

The role of the soft component if further illustrated by fig. |12| where we show the integrated 
gluon density (^Sj) and its soft and hard components Ggoftix^Q"^) and GhardiXjQ"^), respec- 
tively. The soft contribution Gsoft{x, Q"^) is a dominant feature of the integrated gluon density 
Gd{x, Q"^) for <1 GeV^. It builds up rapidly with Q"^ and receives the major contribution 
from the region ~ 0.3-i-0.5 GeV^. Our Ansatz for J^soft{x, k^) is such that it starts decreas- 
ing already at ~ 0.2 GeV^ and vanishes rapidly beyond > k^^j^, see figs. 10,11. Still the 
residual rise of the soft gluon density beyond ~ 0.5 GeV^ is substantial: Gsoft{x, Q"^) rises 
by about the factor two before it flattens at large Q^. We emphasize that GsoftiQ^) being 



20 




Figure 10: D-GRV differential gluon structure function jF(x, k^) as a function of h? at sev- 
eral values of X. Dashed and dotted lines represent the soft and hard components; the total 
unintegrated gluon density is shown by the solid line 

finite at large is quite natural — a decrease of GgoftiQ"^) at large Q"^ only is possible if 
^softi.Q'^) becomes negative valued at large Q^, which does not seem to be a viable option. 

At moderately small x = 10~^ the scaling violations are still weak and the soft contribution 
Gsoft{x,Q^) remains a substantial part, about one half, of integrated GSF Gr){x,Q^) at all 
Q^. At very small x < 10~^ the scaling violations in the gluon structure function are strong 
and GhardiXjQ"^) > Gsoft{x,Q'^) starting from Q"^ ~ 1-2 GeV^. 

6.3 Soft /hard decomposition of the proton structure function ^2(0:, Q^) 

Eqs. (|15|), ( p!6D define the soft/hard decomposition of the proton structure function. In fig. |T3| 
we show F2p'^'^{x, Q^) and F2p^^{x, Q'^) as a function of for the two representative values of 
X. Notice how significance of soft component as a function of rises from fully differential 
T{x^Q'^) to integrated Gd{x,Q'^) to doubly integrated F2p^^{x,Q'^). At a moderately small 
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Figure 11: The same as in Fig. but overlaid onto one graph for illustration of the x- 
dependence of J-'{x, k^). The dashed lines shows the soft component J-'soft{x, k^) and its slight 
variation with x due to the finite-x factor (1 — x)^ . 



a; ~ 10 ^ the soft contribution is a dominant part of F2p{x,Q'^), although the rapidly rising 
hard component F2p^'^{x, Q^) gradually takes over at smaller x. 

Notice that not only does F2p^*{x, Q^) not vanish at large Q^, but also it rises slowly with 
Q"^ as 

T?soft( 2^^soft{Q'^) , 1 

Again, the decrease of F2^°'^*(x, Q^) with would only be possible at the expense of unphysical 
negative valued GgoftiQ^) at large Q^. 

7 DGSF in the x-space: effective intercepts and hard- 
to-soft diffusion 

It is instructive to look at the change of the x-dependence from the differential gluon struc- 
ture function J^{x,Q^) to integrated gluon structure function Gd{x,Q'^) to proton structure 
function F2p{x,Q^). It is customary to parameterize the x dependence of various structure 
functions by the effective intercept. For instance, for the effective intercept r^ff for differential 
gluon structure function is defined by the parameterization 

T{x, K^) oc f - j . (47) 
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Figure 12: The same as Fig. 11 hut for integrated gluon structure function Go{x, Q"^) as given 
by the D-GRV parameterization of the differential gluon structure function jF(x, k"^), for the 
discussion see Section 6.2. 



One can define the related intercepts Thard for the hard component J-'hard{x,Q'^). Notice, that 
in our Ansatz Tgoft = 0. 

The power law ( ^Tj) is only a crude approximation to the actual x dependence of DGSF 
and the effective intercept T(.ff will evidently depend on the range of fitted x. To be more 
definitive, for the purposes of the present discussion we define the effective intercept as 

log(xi/a;2) 

taking X2 = 10~^ and xi = 10^^. The effective intercept ThardiQ"^) is defined by (^81) in terms 

of J^hard{x,Q^). 

One can define the related intercepts Ag//, Xhard for the integrated gluon structure function 
G'd(x,Q2): 

G'b(x,Q2)oc - . (49) 



.X, 

In the case of F2p{x, Q^) we define the intercept A((5^) in terms of the variable x defined 



as 



where My is the mass of the ground state vector meson in the considered flavor channel. Such 
a replacement allows one to treat on equal footing <1 GeV^, where the formally defined 
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Figure 13: The soft-hard decomposition of n-factorization results for the proton structure func- 
tion F2p{x,Q'^) evaluated with the D-GRV parameterization of the differential gluon structure 
function T{x^ k^). 



Bjorken variable can no longer be interpreted as a lightcone momentum carried by charged 
partons. For the purposes of the direct comparison with r(Q^), A(Q^) and in order to avoid 
biases caused by the valence structure function, here we focus on intercepts /S.eff)^hard for 

'2p 



the sea component of the proton structure function F^^'^^x^Q'^) 



F4-(x,g2)cxy . (51) 

The results for the effective intercepts are shown in figs. [T^, [l^ and |l^. 

In our simplified hard-to-soft extrapolation of J-'hard{x, Q^) we attribute to J-'hard{x, Q^) 
at < Q2 ^j^g 

same x-dependence as at = Ql modulo to slight modifications for the 
x-dependence of This gives the cusp in Thard{Q^) at = Q^, i.e., the first derivative of 
ThardiQ"^) is discontinuous at = Ql- 

A comparison of fig. |TT| with fig. and further with fig. |T3| shows clearly that only in 



DGSF J-'{x, Q"^) the effect of the soft component is concentrated at small Q^, in integrated 
Gd{x,Q^) and especially in the proton structure function F2p{x,Q^) the impact of the soft 
component extends to much larger Q^. The larger the soft contribution, the stronger is the 
reduction of Tg// from r^ard and so forth, the pattern which is evident from fig. 14 a to 14 b to 
P^, see also figs. |15| and [Tl 



The change of effective intercepts from differential Q^) to integrated Gd{x,Q^) is 
straightforward, the principal effect is that XhardiQ"^) < ThardiQ"^) and XeffiQ^) < ^effiQ^) 
which reflects the growing importance of soft component in Gd{x, Q^). The change of effective 
intercepts from JF(x, Q^) and Gd{x, Q^) to F2p{x, Q^) is less trivial and exhibits two dramatic 
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Figure 14: Effective intercepts for total, and hard components of, (a) the differential gluon 
structure function jF(x, Q^); (h) integrated gluon structure function Goix^Q"^) and (c) pro- 
ton structure function F2p{x,Q'^) evaluated with the D-GRV parameterization of the differ- 
ential gluon structure function J^{x,h?). In the box (d) we compare the effective intercepts 
'Teff{Q'^)Aeff{Q'^) and AeffiQ'^) for J^{x,Q'^), Gd{x,Q'^) and F2p{x,Q'^), respectively. 

consequences of the hard-to-soft and soft-to-hard diffusion. If the standard DGLAP contribu- 
tion (^) were all, then the change from the intercept A((5^) for integrated gluon density to 
the intercept A(Q^) for the proton structure function F2p{x,Q'^) would have been similar to 
the change from t(Q^) to A(Q^), i.e., the effective intercept Ae//(Q^) would have been close to 
zero for <1 GeV^. However, by virtue of the hard-to-soft diffusion phenomenon inherent 
to the K-factorization, F2p{x, Q"^) receives a contribution from gluons with k,"^ > Q"^, which 
enhances substantially A^ardiQ'^) and A^ffi^Q^). The net result is that at small to moderately 
large we find A^ardiQ^) > ^hardiQ"^) and A^jf^Q^) > \eff{Q^)- As we emphasized above 
in sections 5.3, the rise of real photoabsorption cross section is precisely of the same origin. 

The second effect is a dramatic flattening of effective hard intercept, AhardiQ'^) , over the 
whole range of Q^. For all three DGLAP inputs AhardiQ^) flattens at approximately the same 
Aw ~ 0.4. 

The whole set of flgs. P^|-[T^ also shows that the systematics of intercepts in the hard 
region of > Ql is nearly identical for all the three DGLAP inputs. In the soft region 
we have a slight inequality Thard{i^'^)\D_MRS ^ 'Thard{i^'^)\£)^GRV^ which can be readily at- 
tributed to a slight inequality Ql{MRS) > Ql{GRV). In the case of CTEQ4L(v.4.6) in- 
put the value of QliCTEQ) is substantially larger then Ql{MRS),Ql{GRV). In the range 
Ql{MRS),Ql{GRV) < < Ql{GTEQ) the effective intercept t^^k^) rises steeply with 
K^. This explains a why in the soft region Thard{i^'^)\cTEO signiflcantly larger than for the 
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Figure 15: Effective intercepts for total, and hard components of, (a)the differential gluon 
structure function jF(x, Q^); (h) integrated gluon structure function Gd{x.,Q'^) and (c) pro- 
ton structure function F2p{x, Q^) evaluated with the D-MRS parameterization of the differ- 
ential gluon structure function J^{x,h?). In the box (d) we compare the effective intercepts 
Teff{Q'^)Aeff{Q'^) and AeffiQ'^) for J^{x,Q'^), Gd{x,Q'^) and F2p{x,Q'^), respectively. 



D-GRV and D-MRS parameterizations. The difference among intercepts for tlie tliree pa- 
rameterizations decreases gradually from differential T{x,h?) to integrated GEtix^Q"^) gluon 
density to the proton structure function F2p{x,Q'^). 

Finally, in fig. Owe compare our results for /S.effiQ'^) with the recent experimental data 
from ZEUS collaboration 
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In the experimental fit the range of x = [xmax-, Xmm] varies from 
point to point, in our evaluation of Ae// from eq. ( ^5|) we mimicked the experimental procedure 
taking X2 = Xmax and xi = Xmin- This explains the somewhat irregular dependence. The 
experimental data include both sea and valence components. At > Ql{GRV) = 0.9 GeV^ 
we included the valence component of the structure function taking the GRV98LO parame- 
terization. For CTEQ4L(v.4.6) and MRS-LO-1998 the values of are substantially larger. 
However, the valence component is a small correction and we took a liberty of evaluating the 
valence contribution F^^^{x,Q^) for Ql{GRV) < < Ql{MRS),Ql{GTEQ). The overall 
agreement with experiment is good. Difference among the three parameterization is marginal 



and can of course be traced back to figs. IT3-|T6 
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Figure 16: Effective intercepts for total, and hard components of, (a) the differential gluon 
structure function jF(x, Q^); (b) integrated gluon structure function GoiXjQ"^) and (c) proton 
structure function F2p{x,Q'^) evaluated with the D-CTEQ parameterization of the differen- 
tial gluon structure function J-'{x,K^). In the box (d) we compare the effective intercepts 
Teff{Q'^)Aeff{Q'^) and AeffiQ'^) for J^{x,Q'^), Gd{x,Q'^) and F2p{x,Q'^), respectively. 

8 How the gluon densities of a^- factorization differ from 
DGLAP gluon densities 

It is instructive also to compare our results for integrated GSF (^) with the conventional 
DGLAP fit Gpt^x^Q"^). In fig. |1^ we present such a comparison between our integrated D- 
GRV distribution (the solid curves) and the GRV98LO distribution (the dashed curves). As 
was anticipated in section 3.2, at very large the two gluon distributions converge. We also 
anticipated that at small x and moderate the DGLAP gluon structure functions Gpt{x, Q^) 
are substantially larger than the result of integration of DGSF, see eq. (^). At a; = 10~^ they 
differ by as much as the factor two-three over a broad range of < 100 GeV^. The difference 
between integrated DGSF and the DGLAP fit decreases gradually at large x, and is only 
marginal ai x = 10~^. 

Recall the substantial divergence of the GRV, MRS and CTEQ gluons structure functions 
of DGLAP approximation Gpt{x,Q'^) at small and moderate Q^. Contrary to that, the k- 
factorization D-GRV, D-CTEQ and D-MRS gluon structure functions Gd{x,Q^) are nearly 
identical. We demonstrate this property in fig. ^ where we show integrated GoiXjQ"^) and 
their DGLAP counterparts Gpt{x, Q^) for the three parameterizations at two typical values of 
X. Because of an essentially unified treatment of the region of < Ql and strong constraint 
on DGSF in this region from the experimental data at small Q^, such a convergence of D-GRV, 
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Figure 17: Effective intercepts A(Q^) of the proton structure function F2p{x, Q"^) in the HERA 
domain evaluated for the D-GRV, D-MRS and D-CTEQ parameterizations for the differential 
gluon structure function jF(x, k^); the experimental data points are from ZEUS 



D-CTEQ and D-MRS DGSF's is not unexpected. 

One can also compare the effective intercepts for our integrated GSF Gd{x,Q^) with 



those obtained from DGLAP gluon distributions Gpt{x,Q'^). Fig.EOl shows large scattering of 



A?/)(Q2) from one DGLAP input to another. At the same time, this divergence of differ- 
ent DGLAP input parameterizations is washed out to a large extent in the /^-factorization 
description of physical observables (see also p!7|) . 

9 How different observables probe the DGSF jF(x, Q^) 

The issue we address in this section is how different observables map the k,^ dependence of 
J-'{xg, K^). We expand on the qualitative discussion in section 3.2 and corroborate it with nu- 
merical analysis following the discussion in pOf. . We start with the two closely related quan- 
tities — longitudinal structure function Fl{x,Q'^) and scaling violations dF2{x,Q^)/d\ogQ^ 
— and proceed to F2p{x, Q"^) and the charm structure function of the proton F^pl^x, Q^). This 
mapping is best studied if in (^) and (0) we integrate first over k and z. In order to focus 
on the dependence we prefer presenting different observables in terms of T{2x, k?) and 
G'd(2x, k^) 

F,{x,Q') = ^^E4 / ^eF)(Q^^.^)^(2x,K2), (52) 
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Figure 18: Comparison of our results for integrated gluon density Gd{x,Q^) evaluated with 
the D-GRV parameterization of the differential gluon structure function !F{x,k'^) with the 
GRV98L0 DGLAP input parameterization Gpt{x, Q^). 
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Figure 19: A comparison of the divergence of GRV98L0, CTEQ4L(v.4.6) and MRS-LO-1998 
gluon structure functions Gpt{x, Q^) in the left box with the divergence of our integrated gluon 
structure functions Gd{x, Q^) evaluated for the D-GRV, D-CTEQ and D-MRS parameteriza- 
tions for differential gluon structure function T{x^ Q'^) at two typical values of x 
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Figure 20: A comparison of the intercept ^^//{Q'^) of the x-dependence of the GRV98L0, 
CTEQ4L(v.4-6) and MRS-LO-1998 gluon structure functions Gpt{x^Q'^) with their counter- 
part XeffiQ"^) for integrated GoiXiQ^) evaluated with D-GRV, D-GTEQ and D-MRS param- 
eterizations for differential gluon structure function T{x^ Q"^). 



dF2{x,Q^) _ as{Q^ 
dlogQ'^ Stt 



Qi^^\Q\H?)J^{2x,H?) 



(53) 



In the numerical calculation of Fl{x, Q"^) starting from eq. (|r^) we have Xg and as the 
two running arguments of T{xg,H?). As discussed above, the mean value of Xg is close to 
2a;, but the exact relationship depends on h?. The k, z integration amounts to averaging of 
T[xg^ K^) over certain range of Xg. The result of this averaging is for the most part controlled 
by the effective intercept rejj(«;^): 



{Tixg.n})) = ( J^(2x, 



2x 



X n 



(54) 



Because the derivative of Teff{n, ) changes rapidly around k = Q^, the rescaling factor r(K ) 
also has a rapid variation of the derivative at = Q^, which in the due turn generates 

the rapid change of derivatives of 9£^('*(Q^, /t^) around = Q'^. As far as the mapping of 
differential J-'{2x, h<?) is concerned, this is an entirely marginal effect. However, if we look 
at the mapping of integrated gluon structure function Gd{x, Q"^), which is derived from (|52D , 
33f) by integration by parts: 



3n 



dK,^ dQ^P 



{Q 



2 



d log 



Gd{2x,h? 



(55) 



dF2{x,Q^ 
d log 



37r 



E4 



kP' d log 



-Gd{2x,h? 



(56) 
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Figure 21: The weight function for mapping of the differential gluon structure function 
jF(x, K^) as a function of for several values of Q"^ . We show separately the results for light 
flavours, u, d, and charm. 

then the weight functions (9G2'^/^(Q^, K^)/91og/t^ will exhibit a slightly irregular behaviour 
around k;^ = Q"^. Evidently, such an irregularity appears in any region of fast variation of 
Tf.ff{K^)] in our simplified model it is somewhat amplified by the cusp-like dependence of 

Finally, starting from ( |5B| ) one obtains a useful representation for how the proton structure 
function F2p(x, Q"^) maps the integrated gluon structure function: 

F ix O') - V e' ( 9_^i!lSllJ^c i2x 



= / ^wi''~\Q',^')asi>^')Gni2x,^') (57) 

In figs. and ^ we show the weight functions and B2. Evidently, for light flavours 
and very large they can be approximated by step-functions 

Q^li\Q',k')-^9{CL,2Q'-i^'), (58) 

where the scale factors Cl ~ | and C2 ~ 2 can be readily read from figures, for the related 
discussion see [^. Note that the value C2 ~ 2 corresponds to 6*2 ~ 8 introduced in Section 3.2. 



Recall that the development of the plateau-like behaviour of and G2 which extends to ~ 
signals the onset of the leading logQ^ approximation. For large in the approximation 
(|58D the integration can be carried out explicitly and Fl{x, Q"^) oc Gd{2x, ClQ^)- Similarly, 
dF2{x,Q^)/dlogQ^ oc GDi2x,C2Q^), cf. eq. (H). 
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Figure 22: T/ie weight function Q2 for mapping of the differential gluons structure function 
jF(x, K^) as a function of h? for several values of . We show separately the results for light 
flavours, u, d, and charm. 



Still better idea on how Fl and scaling violations map the integrated GSF is given by 
figs. 0,0, where we show results for —dQ{ff)j^/d\ogH? and W2^\ The first quantity is 
sharply peaked at ClQ^. The second quantity visibly develops a plateau at large Q^. As 
can be easily seen, scaling violations do receive a substantial contribution from the beyond- 
DGLAP region of h?>Q'^. 

Because of the heavy mass, the case of the charm structure function F2p{x, Q"^) is somewhat 
special. Figs. [0 and ^ show weak sensitivity of F2p{x, Q^) to a soft component of J-'{xbj, k,"^) 
which has an obvious origin: long- wavelength soft gluons with k < rric decouple from the color 
neutral cc Fock state of the photon which has a small transverse size < Our results for 
F2p{x, Q^) are shown in fig. ^ the agreement with the recent precision experimental data 
from ZEUS is good. 



10 Summary and outlook 

We present the first parameterization of differential gluon structure function jF(x, Q"^) of the 
proton inherent to the K-factorization approach to small-x DIS. The form of the parameteri- 
zation is driven by color gauge invariance constraints for soft Q^, early ideas from color dipole 
phenomenology on the necessity of nonperturbative soft mechanism for interaction of large 
color dipoles and by matching to the derivative of familiar DGLAP fits Gpt{x, Q"^). The latter 
condition is not imperative, though, and can be relaxed; in this exploratory study we simply 
wanted to take advantage of the insight on Gpt{x, Q"^) from early DGLAP approximation stud- 
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Figure 23: The same as Fig. |^ but for mapping of the integrated gluon structure function 
Gd{x, k^) as a function of for several values of . We show separately the results for light 
flavours and charm. 




Figure 24: The weight function W2 for mapping of the integrated gluons structure function 
Gd{x, k^) as a function of for several values of . We show separately the results for light 
flavours and charm 

ies on scaling violations. The parameters of Q'^) have been tuned to the experimental 
data on in the low-x (x < 0.01) domain and throughout the entire region as well as on 
real photoabsorption cross section cr/jj- Differential gluon structure function jF(x, Q^) is the 
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Figure 25: A comparison of the experimental data from ZEUS / p^/ on the charm structure 
function of the proton with n-factorization results for F2'^{x, Q^) based on the D-GRV param- 
eterization of the differential gluon structure function J-'{x, Q"^). 

principal input for pQCD calculation of many diffractive processes and we anticipate that the 
consistent use of our parameterizations shall reduce the uncertainties of calculations of cross 
section of such processes as diffractive DIS into vector mesons and continuum. 

Our results allow to address several interesting issues. First, our Ansatze for J-'{x, Q^) have 
been so constructed as to ensure the convergence of Gd{x, Q^) — the integral of J-'{x, Q"^) — 
to the corresponding large DGLAP input Gpt{x,Q^). We notice that both gluon distribu- 
tions provide a comparable description of the same set of the experimental data on the proton 
structure function, the only difference being that in the K-factorization we lift the DGLAP lim- 
itation on the transverse phase space of quarks and antiquarks. We find very slow convergence 
of, and numerically very large difference between, the /t-factorization distribution GoiXjQ"^) 
and the DGLAP fit Gpt{x,Q'^). As anticipated, the divergence of the two distributions is 
especially large at small-x and persists even in the hard region up to ~ 10 -i- 100 GeV^ at 
X = 10~^. We interpret this divergence as a signal of breaking of the DGLAP approximation 
which arguably gets poorer at smaller x. The second finding is a numerically very strong 
impact of soft gluons on the integrated gluons structure function Gd^x^Q"^) and the proton 
structure function F2p{x,Q'^). It is not unexpected in view of the early work on color dipole 
phenomenology of small-x DIS, but the evaluation of the soft component of integrated gluon 
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structure function is reported here for the first time. In conjunction with the strong departure 
of the K-factorization distribution Ge>{x,Q'^) from the DGLAP fit Gpt{x,Q'^) it serves as an 
warning against unwarranted apphcation of DGLAP evolution at in the range of several 
GeV^. 

The phenomenologically most interesting finding is the anatomy of the rising component 
of the proton structure function from the Regge theory point of view. We notice that effective 
intercepts ThardiQ"^) and \hard for hard components of the differential and integrated gluon 
distributions are lively functions of which vary quite rapidly with from ?s 0.1 at small 

to ~ 0.6 at ~ 10^ GeV^. In the Regge theory language this evidently implies that hard 
component of neither jF(x, Q^) nor Gd^x^Q"^) is dominated by a single Regge pole exchange 
and a contribution from several hard Regge poles with broad spacing of intercepts is called 
upon. However, an approximately fiat dependence of /^hardiQ"^) shows that the hard 
component of the proton structure function can well be approximated by a single Regge pole 
with intercept ^hard ~ 0.4. Such a scenario in which a contribution of subleading BFKL- 
Regge poles to F2p{x, Q"^) is suppressed dynamically because of the nodal properties of gluon 
distributions for subleading BFKL- Regge poles has been encountered earlier in the color dipole 
BFKL approach The intercept A^ardiQ'^) found in the present analysis is remarkably close 
the intercept of the leading BFKL- Regge pole Ajp = 0.4 found in the color dipole approach in 



1994 p6| , |l^, |2T[. For the related two-pomeron phenomenology of DIS see also [0. From the 
point of view of K-factorization, the hard-to-soft diffusion is a unique mechanism by which an 
approximate constancy of AhardiQ"^) derives from a very rapidly changing ThardiQ"^)- Fourth, 
the same hard-to-soft diffusion provides a mechanism for the rise of the real photoabsorption 
cross section a''^ in a model with the manifestly energy independent soft cross section. We 
emphasize that the hard-to-soft diffusion is a generic phenomenon and we do not see any 
possibility for decoupling of hard contribution from photoabsorption at = 0. 

We restricted ourselves to a purely phenomenological determination of differential gluon 
distributions from the experimental data on F2p{x,Q'^) which is sufficient for major appli- 
cations of the At-factorization technique. Whether the so-determined hard components of 
Q^) and Gd^x^Q"^) do satisfy the dynamical evolution equations and what is the onset 
of DGLAP regime will be addressed elsewhere. 

IPI wishes to thank Prof. J.Speth for his hospitality at Forschungszentrum Jiilich. This 
work was partly supported by the grant INTAS 97-30494. 
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